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ABSTRACT 

We present numerical studies of S-dimcnsional electron magnetohydrodynamic (EMHD) turbulence. 
We investigate cascade timescale and anisotropy of freely decaying strong EMHD turbulence with 
zero electron skin depth. Cascade time scales with k~'^/'^. Our numerical results clearly show scale- 

1/3 

dependent anisotropy. We discuss that the observed anisotropy is consistent with fc|| ex k^ , where 



±_ are wave numbers parallel and perpendicular to (local) mean magnetic field, respectively. 
Subject headings: MHD — turbulence — acceleration of particles 



fell and k 



1. INTRODUCTION 

Electron magnetohydrodynamics (EMHD) deals with 
MHD phenomena occurring at scales smaller than those 
of conventional MHD (see Kingsep, Chukbar, & Yankov 
1990). On scales below the ion inertial length di — c/uipi, 
where c is the speed of light and uopi is the ion plasma fre- 
quency, we can assume that the ions create only smooth 
motionless background and fast electron flows carry all 
the current, so that 

Ve = -— - --^V X B. (1) 

ne 47me 

When we ignore the ambipolar diffusion, the Ohm's law 
becomes E = — Vg x B/c + 77J, where 77 is the plasma 
resistivity. Inserting this into Faraday's induction law 
{dB/dt = -cV X E), we obtain the EMHD equation 

^ " '^e^ "" f^^ X B) X B] + T^V^B. (2) 

Note that, in this paper, we only consider the zero (nor- 
malized) electron inertial length case: d^ = c/{LOp!,L) — > 
0, where ojpe is the electron plasma frequency and L is 
the typical size of the system. 

EMHD turbulence (also known as whistler turbulence) 
can be viewed as a generalized Alfvenic turbulence on 
scales smaller than the proton Larnror radius (Quataert 
& Gruzinov 1999). EMHD turbulence plays important 
roles in neutron stars and accretion disks. Goldreich & 
Reisenegger (1992) discussed the properties of EMHD 
turbulence and argued that EMHD turbulence can en- 
hance ohmic dissipation of magnetic field in isolated neu- 
tron stars (see also Cumming, Arras, & Zweibel 2004). 
Quataert & Gruzinov (1999) discussed transition from 
conventional Alfvenic MHD cascade to EMHD cascade in 
advection dominated accretion flows (ADAF). Other ap- 
plications of the EMHD model include quasi-coUisionless 
magnetic reconnection in laboratory and space plas- 
mas (Bulanov, Pegoraro, & Sakharov 1992; Biskamp, 
Schwarz, & Drake 1995; Avinash et al. 1998) and plasma 
opening switches and Z pinches (see Kingsep et al. 1990). 

Earlier researchers convincingly showed that energy 
spectrum of EMHD turbulence is steeper than Kol- 
mogorov's k^^^^ spectrum. Using two-dimensional (2D) 



numerical simulations, Biskamp, Schwarz, & Drake 
(1996) found that energy spectrum follows E{k) oc fc"'' 
with n = 2.25 ± 0.1. They showed that the following 
Kolmogorov-type argument can explain this spectrum. 
Suppose that the eddy interaction time for eddies of size 
/ is the usual eddy turnover time tcas,i ~ l/vi- Since 
V ex V X B (equation (^)), this becomes tcas,i oc P/bi. 
Combining this with the constancy of spectral energy 
cascade rate {bf /tcas.i—con.sta.nt), one obtains E{k) ex 
k^'^/^. Note that E{k) and k are related by kE{k) ~ bf. 
Biskamp et al. (1999) also obtained a similar result for 
3D. Ng et al. (2003) confirmed the 2D result of Biskamp 
et al. (1996). Therefore, earlier works strongly support 
that 

E{k) ex k-"^/^ (as 4 -> 0) (3) 

for both 2D and 3D. Recently, Galtier & Bhattacharjee 
(2003) studied energy spectrum of weak whistler turbu- 
lence. 

Dastgeer et al. (2000), Ng et al. (2003), and Dastgeer & 
Zank (2003) numerically studied anisotropy of 2D EMHD 
turbulence. No 3D study has been reported yet. In this 
paper, we study anisotropy of 3D EMHD turbulence. 

Method. — We have calculated the time evolution of 
decaying 3D incompressible electron magnetohydrody- 
namic turbulence. We have adopted a pseudospectral 
code to solve the normalized incompressible EMHD equa- 
tion in a periodic box of size 2tt: 

dB 

'dt 

where magnetic field, time, and length are normalized by 
a mean field Bq, the whistler time i^, = L^(wpe/c)^/r2e 
{fle= electron gyro frequency), and a characteristic 
length scale L (see, for example, Galtier & Bhattachar- 
jee 2003). The resistivity rj' in equation Q is dimen- 
sionless. The dispersion relation of a whistler waves in 
this normalized units is lu — kknBo, where ku is the wave 
number parallel to the (local) mean magnetic field. The 
magnetic field consists of the uniform background field 
and a fiuctuating field: B — Bq + b. The strength of the 
uniform background field, Bq, is set to 1. We use 288'^ 
collocation points. At t — 0, the random magnetic field 



= -V X [(V X B) X B] + Ty'V^B, 



(4) 
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- Energy spectra, timescale and anisotropy. (a) Energy spectra for t=0 (dotted), 0.33 (dashed), and 0.54 (solid). The spectra 
33 are compatible with E{k) oc k~'^'^. (b) Cascade time scales with teas oc A;~^'^. (c) Anisotropy calculated from equation 171 . 



is restricted to the range 2 < fc < 4 in wavevector (k) 
space. The ampHtudes of the random magnetic field at 
i = is '^ 1. Hyperdiffusivity is used for the diffusion 
terms. The power of hyperdiffusivity is set to 3, so that 
the dissipation term in the above equation is replaced 
with 773(V^)^B, where 773 is approximately 2 x 10~^°. 

2. RESULTS 

Spectrum. — Figure ^ shows energy spectra. At t = 
only large scale (i.e. small fc) Fourier modes are excited. 
The dotted curve in Figure^shows the initial spectrum. 
As the turbulence decays, the initial energy cascades 
down to small scales and, as a result, small scale (i.e. 
large k) modes are excited. When the energy reaches 
the dissipation scale at fc > 70, the energy spectrum de- 
creases without changing its slope (the dashed and the 
solid curves). The slope at this stage is around 2.20, 
which is very close to the predicted spectrum: 



E{k) oc fc-^/3. 



(5) 



However, we also note that our result is compatible with 
a fc~^ spectrum. 

Timescale. — The energy spectra confirm that the 
Kolmogorov-type energy cascade model with teas '^ l/vi 
works fine for EMHD turbulence. Here we show that 
timescale of motion also supports the energy cascade 
model. 

Symbolically, we can rewrite the EMHD equation in 
wave-vector space as follows: bk — N^, where N^ rep- 
resents the nonlinear term. We have ignored the dissipa- 
tion term. Naively speaking, we might obtain the time 
scale by dividing |bk| by |N^|. However, to remove the 
effect of large scale translational motions, we restrict the 
evaluation of the nonlinear terms to contributions coming 
from the interactions between the mode at fc and other 
modes within the range of fc/2 and 2fc (see Cho, Lazar- 
ian, & Vishniac (2002a) for details). We show the result 
in Figure 1(b). Our result seems to follow the timescale 



teas - I/Vl - ;"/' 



-4/3 



(6) 



Anisotropy. — Since B ■ VB ^ Bq < fc|| > 6, we can 
obtain the average parallel wavenumber < fcii > from 
< fc|| >- (|B ■ VBp/Sg/62)i/2^ where b is the strength 
of the random magnetic field. The quantity 1/ < fcii > is 



roughly the average parallel size of eddies at the energy 
injection scale. We can easily extend this method to 
smaller scale eddies. Consider eddies of perpendicular 
size I. Since it is reasonable to assume that the eddies 
are stretched along the local mean magnetic field, the 
quantity Bi-Vb/, not Bo-Vb/, is proportional to -Bifcy 5/, 
where B^ is the local mean field, b; the fluctuating field 
at scale I, and fcn the wave number parallel to the local 
mean magnetic field (Cho, Lazarian, & Vishniac 2002b). 
Therefore, we have fc||(0 - (|Bl ■ Vhif / Bl/bfY/'^ . It is 
convenient to do the calculation in Fourier space: 



fc||(fc_L) 



E 



fe<|k'|<fe-|-l 



IBi-VbH^, 



1/2 



^L Z^A;<|k'|<fcH 



1 Ib|^. 



(7) 



where fcj_ ('--^ fc ^^ l/l when anisotropy is present) is the 
wave number perpendicular to the local mean magnetic 
field and Bl is the local mean field obtained by eliminat- 
ing modes whose perpendicular wavenumber is greater 
than fc/2 (~ 1/(2Z)). The fluctuating field b; is obtained 
by eliminating modes whose perpendicular wavenumber 
is less than fc/2. We plot the result in Figure 1(c). The 
result does not show a well-defined power law between fc|| 
and fcj_. However, the result is roughly compatible with 



fcii 



.1/3 



(8) 



•|| --i 

We believe this is the true anisotropy of EMHD turbu- 
lence. 

We may use the second order structure function to il- 
lustrate scale-dependent anisotropy of eddy structures 
(see Cho & Vishniac 2000). We calculate the second or- 
der structure function in a local frame, which is aligned 
with local mean magnetic field Bl: 

SF2(r|| , r^) =< |B(x + r) - B(x)P >avg. over X, (9) 

where r — r||f|| -I- rxf± and f || and fx are unit vectors 
parallel and perpendicular to the local mean field Bl, 
respectively. See Cho et al. (2002a) and Cho & Vishniac 
(2000) for the detailed discussion of the local frame. We 
show the resulting structure function in Figure 2. 

We may obtain anisotropy of eddy structures by an- 
alyzing the contour shapes. However, we believe that 
this method will not give the true anisotropy for EMHD. 
The reason is as follows. Roughly speaking, when en- 
ergy spectrum is steeper than fc^^ we should not use the 
method. When energy spectrum is steeper than fc~^. 
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Fig. 2. — Contour plot of the second structure function 
SF2{ri\,r^) at i = 0.44, where rn (rj^) is distance along (perpen- 
dicular to) the local mean field, shows existence of scale-dependent 
anisotropy: smaller eddies are more elongated. Although this plot 
is useful to illustrate scale-dependent anisotropy, we do not use this 
plot to study anisotropy (see text). 



small scale fluctuations are so small that the calculation 
of the second-order structure function will be dominated 
by the smooth change of the large scale field and the 
second-order structure function will show r^-scaling re- 
gardless of the true scaling: < \B{x) — B{x + r)p >~< 
\Bl{x) — Bl{x + r)\'^ >oc r^ when r < L. In the perpen- 
dicular direction, E{k±) oc fcj^ and, therefore, it is fine 
to use the structure function. However, we should not use 
the method for the parallel direction when anisotropy is 

2/3 

stronger than fcii oc fc^ . Suppose anisotropy scales with 
k 

k 



oc fc" . Then the 3D energy spectrum is £^3/3 (fcn ,kj_) on 



II 
-10/3-Q 



.g(fc||/fc"), where g is a function that describes 



distribution of energy along the fcii direction in Fourier 
space (see Goldreich & Sridhar 1995; Cho et al. 2002a). 
The ID energy spectrum for the parallel direction be- 
comes £'(fc||) oc J kj^d(j)dk±E3£i{k\\,ki_) oc k"., 
(Cho et al. 2002a). Therefore, when a < 2/3, £;(fc||) 
is steeper than fc7 and, thus, we cannot use the struc- 
ture function to study anisotropy. In other words, even if 
a < 2/3, the structure function method will give a wrong 
result that a « 2/3. 
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Fig. 3. — Typical eddy shape and bi/Bj^ ratio. The degree of 
bending of magnetic field lines depends on the bi/Bi^ ratio. When 
wo assume that lj_ and lu represent typical perpendicular and par- 
allel size of eddies at the scale I (~ lj_), we obtain l±/lii ~ 6; /-Bo 
and, thus, Zn '^ I / , we use 6; '^ P'^ and B^ ~ -Bo- 



along the magnetic field line at the speed of kBo. When 
this wave packet collides with opposite-traveling wave 
packets (of similar size), the change of energy per colli- 
sion is AE - {dE/dt)At - {k'^bf){li^/{kBo)), where we 
used equation Q) to estimate dE/dt. Therefore, 



AE/E ~ (fc6/)/(A:||So) - t^ /teddy 



(10) 



If AE/E < 1, then many collisions are required to make 
AE/E ^ 1. These collisions are incoherent. Then 
{teddy /twY collisions are necessary to complete the cas- 
cade. This means that the energy cascade time by colli- 
sion is 



'^cas^coll ^^ y^eddy / '^w ) ^^ 



(11) 



If teas, coll ~ teddy, this equation implies t^ ~ teddy 
Therefore, when teas ^ teddy, the turbulence is strong 
(i.e. AE/E '^ 1) ^ and there is a balance between 
teddy ~ ^/{k'^h) and tw ~ l/(A:A:||-Bo)- Our results 
clearly show that the critical balance factor ^ = t^ /teddy 
is constant (note that bi ~ P^^ and ku oc fc^/'^). Us- 



ing bl 



k •^/'^, we obtain fcn 



oc 



.1/3 



The 



tainty relation is another way to express this balance. 
The usual uncertainty relation between time and energy 
(^frequency) is At oc l/Aw. Since Ai ~ teas ^ k~'^/'^ 
and Auj ^ B^kkn , we obtain the same anisotropy scaling. 
Second, a simple geometric argument in Fig.Olsuggests 



We obtained k 



3. DISCUSSION 
X kl in the previous section. See the fcii oc fcV^ anisotropy. Suppose that Fig. 01 shows a 



equation ^. Here we discuss why this is reasonable. 

First, the fact that the Kolmogorov-type argument 
works fine for the energy spectrum implies the k\\ oc 

1/3 

fc_/ anisotropy. As we mentioned earlier, Biskamp et 
al. (1996) assumed that the cascade time {teas) is the 
eddy turnover time {teddy ^ ^/vi ^ I'^/h) and obtained 
E{k) oc k~'^'^, which agrees well with numerical simula- 
tions. See also our calculation of timescale in Figure 1 (b) . 
We can show that, when teas ^ teddy, there is balance be- 
tween the eddy turnover time teddy and the the whistler 
time i„, (=the linear wave period=^ l/(fcA;||-Bo)). Sup- 
pose that we have a wave packet whose parallel size is 



typical shape of magnetic field lines at scale /. Then, we 
can assume that l±/l\\ ~ hi/B^, where 6; is the fluctu- 
ating field at the scale I and Bl is the local mean field 
{Bl ~ Bq). If we assume l±_ and ^y represent typical 
shape of eddies at the scale / (^ ^_l), we have 



kn/k^r^bi/Bo. 



(12) 



fc|| and perpendicular size l± 



when anisotropy is present). 



This wave packet travels 



^ The reverse is also true: if turbulence is strong, then AE/E ^ 
1 by definition and, therefore, energy cascade occur s du ring one 
wave period, which is equal to ieddy (see equation IIUI ). Thus, 
one can show that three expressions, Kolmogorov-type scaling for 
cascade time {teas ~ l/^i), critical balance (t^ddy '^ tw), and strong 
turbulence {AE/E ~ 1), are all equivalent (C. Thompson, private 
communication), as in ordinary MHD turbulence. See Goldreich 
& Sridhar (1995) for the original discussion of the critical balance 
and anisotropy of ordinary MHD turbulence. 
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In this interpretation, smaller eddies are more elongated 
because they have a smaller bi/Bo ratio. Note that the 
relation in equation 1)12(1 is in fact equivalent to the strong 
turbulence condition {AE/E ^ 1; see equation ()10|l '). 
which means energy cascade occurs over one wave period. 
From equation ((l()|l . we can show that teas ~ l/vi ^ l"^ /bi. 
Combining this with the constancy of spectral energy 
cascade rate (&j^/icas=constant), we obtain bi oc k~'^/^ 

and, thus, the k\\ oc k^ anisotropy. 

One peculiar feature of EMHD cascade is that there 
is balance between parallel and perpendicular cascade 
by self-interaction. The resonance relations of whistler 
turbulence, 



ki+k2 = k3 and fci|fc||4l + fc2|fc|| 2I = fcal^ 



l,3l 



(13) 



allow self-interaction. We can show that the parallel cas- 
cade by self interaction is as fast as the perpendicular 



cascade when k 



k 



1/3 



Ng et al. (2003) performed 



a numerical experiment and found that cascade by self- 
interaction is as fast as that by collision with opposite- 
traveling wave packets. Combined together, these facts 

1/3 
imply that anisotropy cannot be stronger than k\\ oc k^ 

in EMHD. We will discuss this in detail in a forthcoming 
paper. 

We have focused on anisotropy of 3D EMHD turbu- 
lence. However, it is fair to mention earlier 2D results. 
Dastgeer et al. (2000) numerically studied 2D EMHD 
turbulence and found numerical evidence that the turbu- 
lence is anisotropic in the presence of a mean field. Dast- 
geer & Zank (2003) showed that the global anisotropy of 
decaying EMHD turbulence scales linearly with the in- 
verse of mean magnetic field strength, \/Bq. Ng et al. 

(2003) found local anisotropy of k\\ oc k/^ from a 2D 
EMHD simulation. We will discuss the 2D versus 3D 
case elsewhere. 



Conventional MHD turbulence can be decomposed into 
cascades of Alfven, fast and slow modes (Goldreich & 
Sridhar 1995; Cho & Lazarian 2002). While fast and 
slow modes get damped at larger scales, Alfvenic modes 
can cascade down to the proton Larmor radius scale. 
These Alfven modes excite whistler modes below the 
scale. Quataert & Gruzinov (1999) discussed transition 
from conventional Alfvenic MHD cascade to EMHD cas- 
cade in ADAFs and argued that anisotropy of EMHD 
turbulence is important for electron heating in ADAFs. 
They showed that, depending on the anisotropy of the re- 
sulting whistler turbulence, the whistler cascade will or 
will not transfer most of its energy to electrons. Whether 
or not whistler turbulence is anisotropic is an important 
question for the physics of accretion disks and gamma- 
ray bursts. The strong anisotropy of EMHD turbulence 
makes heating of protons by EMHD turbulence rather 
difficult in astrophysical plasmas with /3 ~ 1, where [3 
is the gas to magnetic pressure ratio (see Quataert & 
Gruzinov 1999). 

4. SUMMARY 

Using a numerical simulation, we have studied 3D 
EMHD turbulence. We have found that energy cascade 
time scales with k~^^^. We have calculated anisotropy 



and found k 



k 



1/3 



We have discussed why this is 



reasonable. 
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